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Abstract 

We have simulated atmospheric effects such as fog and smoke in laboratory environment 
to simulate depolarisation due to atmospheric effects during a free space optical communi¬ 
cation. This has been used to study noise in two components of quaternary encoding for 
polarization shift keying. Individual components of a Quaternary encoding, such as vertical 
and horizontal as well as 45° and 135°, are tested separately and indicates that the depo¬ 
larization effects are different for these two situation. However, due to a differential method 
used to extract information bits, the protocol shows extremely low bit error rates. The 
information obtained is useful during deployment of a fully functional Quaternary encoded 
PolSK scheme in free space. 


1. Introduction 

Two-binary, one-quaternary (2B1Q) en¬ 
coding is a method of mapping a pair of 
bits (DiBIT) to single encoding of a 4- 
level scheme. Such methods of encoding 
multiple bits to each level of the encoding 
scheme allows increase in transmission den¬ 
sity, although most of the times requires 
hardware with better precision. Several 
schemes for such multiple encoding have 
been proposed and studied earlier III El E] 
and their robustness have been discussed. 
Some of them are optical methods involv¬ 
ing light, wherein different polarization or 
phase states of light are mapped to DiBITs 
PI 1- However, when such methods are 
operated in free space communication pro¬ 
tocols, concern has to be taken about the 
fact that atmospheric phenomena such as 
fog and smoke cause multiple scattering 
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of the light, leading to a depolarization 
H El El El [TO] . In additon, as we show in 
this communication, the depolarization be¬ 
haviour is different for different encodings, 
adding additional difficulty, which is pecu¬ 
liar to M-ary encodings. But our method 
of differential measurement, as shown ear¬ 
lier ra, will allow us to extract informa¬ 
tion with a near zero bit error rate despite 
a significant depolarization. 

The paper is organized as follows. We 
first briefly describe our experimental setup, 
which is same as the one used in mi. We 
then present the results for errors due to 
depolarization, first for vertical and hori¬ 
zontal encoding and then independently for 
45° and 135°. The present study is not a 
complete and proper Quaternary scheme, 
but instead a depolarization study of indi¬ 
vidual components. But this will help us 




understand whether or not the atmospheric 
effects such as fog and smoke affect these 
states differently, and if so, whether that in¬ 
formation can be incorporated into the com¬ 
munication scheme for better detection. 

We also explain the differential method 
of measurement based on State of Polariza¬ 
tion, which allows us a higher tolerance for 
depolarization. Finally we show an analy¬ 
sis of this technique in terms of the Muller 
matrix method, which incorporates the de¬ 
polarization effects. This allows us to repre¬ 
sent the depolarized light as a Stokes Vector 
with partial polarization. The mathemati¬ 
cal analysis for the differential method of 
measurement shows that the final measure 
is a value whose sign can be used to identify 
and map to the information bit. 

2. Experiment 

We consider a specific Quaternary encod¬ 
ing consisting of two pairs of mutually or¬ 
thogonal polarizations - viz Horizontal and 
Vertical, as well as 45°and 135°. These can 
be represented on a standard constellation 
diagram as shown in figure |l(a)[ On the 
left is only the Binary encoding using only 
Horizontal and Vertical polarizations, while 
the right side shows the Quaternary encod¬ 
ing. Angle between them. When finally de¬ 
ployed, we should be able to use this as an 
encoding for pair of bits as 


bit pairs 

Polarization 

00 

vertical 

01 

45° 

10 

135 ° 

11 

horizontal 


Table 1: Quaternary encoding of bit pairs to polar¬ 
ization 



(a) 


Transmitter Atmospheric channel Receiver 



(b) 

Figure 1: (a) Constellation diagrams for Bi¬ 

nary (left) and Quaternary encoding (right) (b) 
Schematic of the experimental setup. See text for 
details. 

The experimental setup is similar to the 
one explained in our earlier paper mi, 
briefly explained here for sake of complete¬ 
ness. It consists of two lasers LI and L2, 
both VCSEL’s operating at 780nm. The 
choice of the wavelength is due to the fact 
that atmospheric attenuation has a clear 
window in this region na. Light from lasers 
LI and L2 is split into two polarized com¬ 
ponents by the polarizing beamsplitter PBS 
and vertical part of light from LI and hor¬ 
izontal part of light from L2 goes into the 
communication channel. A Lab VIEW soft¬ 
ware controls the lasers through the driver 
circuit, such that bit ‘O’ results in a pulse 
from LI and bit ‘1’ results in a pulse from 
L2. The Halfwave plate, indicated by A/2 
in figure, is used to rotate the polariza¬ 
tion whenever required. This is used to 
choose between vertical/horizontal scheme 
or 45°/135°scheme as required. 
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A glass chamber placed in the com¬ 
munication channel is filled with fog or 
smoke to simulate the atmospheric condi¬ 
tions. Smoke is created in the chamber by 
burning household incense powder. Fog is 
achieved by sprinkling water onto a a sam¬ 
ple dry ice kept within. Both smoke and fog 
cause depolarization of the light due to mul¬ 
tiple scattering by smoke particles or water 
droplets. This results in a loss of informa¬ 
tion. 

The amount of smoke or fog is quantified 
by measuring the attenuation of the laser 
after passing through and characterized by 
optical density in dB, as given by [] 


OD(dB) = 10 log 


/ Transmitted Intensity \ 
\ Incident intensity J 


The receiver consists of another polar¬ 
izing beamsplitter, labelled PBS2 and two 
Avalanche Photodiodes (labeled APD1 and 
APD2, both PCD Mini 0020 module from 
SenSL with 20 /rm sensor and pettier 
cooler). The TTL pulses produced by the 
APD’s are recorded via a DAQ card and 
counted using LabVIEW program. The 
counters are synchronized to a stream of 
clock pulses from the transmitter. To char¬ 
acterize the entire system for all possible 
polarization schemes, a random sequence of 
ones and zeros were generated by the com¬ 
puter and transmitted. The transmitted 
and received data are compared to obtain 
true Bit Error Rate (BER). These bit error 
rates are compared with the theoretical es¬ 
timate of the BER for normal OOK scheme. 

The LabVIEW program computes State 
of Polarization, which is defined by us as 


1 An alternative definition of Optical density is in 
terms of Beer - Lambert law, as a = IxiI/Iq, which 
serves the function equally well. In this manuscript, 
we continue to use the definition as in attenuation 


APD2 - APD1 
“ APD2 + APD1 

where APD1 and APD2 indicate the 
counts of photons on the two detectors. In 
the ideal case with no depolarization, a hori¬ 
zontally polarized light will all reach APD1, 
resulting in an SOP of —1, and a verti¬ 
cally polarized light will all reach APD2, 
resulting in an SOP of +1. However, de¬ 
polarization causes some ‘leakage’ to the 
other APD, resulting in an SOP value lesser 
than the range ±1. Since multiple scatter¬ 
ing are random events, each pulse of sim¬ 
ilar polarization results in different values 
of SOP. A distribution of such several SOP 
values obtained from a single run of about 
10000 pulses is shown in figure [2| It clearly 
shows the effect of depolarization, wherein 
the mean values are shifted inward from ±1. 
But since the two distributions are not over¬ 
lapping each other, an unambiguous bit des¬ 
ignation can still be performed, which shows 
a higher tolerance for noise, as shown in our 
earlier work m- We have also shown in 
the earlier work that the measure of SOP 
obtained, is actually due to due to the bal¬ 
listic and snake photons mmm whereas 
those that undergo significant depolariza¬ 
tion cancel out each other. In this com¬ 
munication, we test this method with other 
polarizations, which will form components 
of a Quaternary encoding. 



Figure 2: Statistical Distribution of SOP values 
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3. Results & Discussions 


( 2 ) 


The experiment was performed by trans¬ 
ferring about 1,00,000 bits of random se¬ 
quence of 0’s and l’s, first in Verti¬ 
cal/Horizontal basis (also called a ‘+’ basis) 
and later in 45°/135°(‘X’ basis). Individual 
photon counts from APD1 and APD2 are 
recorded and the State of Polarization, SoP 
was computed as per formula [lj The in¬ 
formation bit was taken to be zero if SoP 
value was negative and bit was considered 
1 if SoP was positive. By transmitting only 
one of the polarizations, we could compute 
the amount of leakage, obtain a histogram 
of the counts for true signal as well as leak¬ 
age (see figure [3]) . 




(a) 


(b) 




(d) 


Figure 3: Histograms of photon counts by APDs. 
(a) is true counts when LI is on and (b) is cor¬ 
responding leakage. Similarly (c) is corresponding 
leakage when only laser L2 is on. 


By fitting the statistical distribution of 
these data to individual Gaussian function, 
we obtain mean photon number for signal 
and noise/leakage as well as width of distri¬ 
bution. Using this in the standard formula 
for Quality factor mm 5 as 


(Jl + cr 0 ’ 

As in case of ON-off keying, I\ and Jo 
stand for mean counts of signal and noise, 
given by centres of Gaussian distribution. 
<7i and do are the width of respective Gaus¬ 
sian distributions. But in this case there 
will be independent Q factors, one each for 
pair polarizations. 


Qi = U Ij 


Oi 


(J i 


where z, j 
45°/135°. 


are 


Qj — 


h ~ h 

Vj + Vi 


( 3 ) 


respectively V, H or 


The overall Q factor would be average of 
the individual Q factors above. The theo¬ 
retical estimate of Bit error rates, as given 
by [TBj , for all possible i 


BER . = erfc(<^V2) 

2 

3.1. The Quality factor 

For the ‘+’ basis transmission, about 
10,000 data bits in random sequence of ze¬ 
ros and ones were transmitted and their in¬ 
dividual SOP values were recorded, as per 
equation [l] Their statistical distribution 
was recorded as a histogram and fit to a 
Gaussian. From the parameters of Gaus¬ 
sian, Q factors were computed and a theo¬ 
retical BER obtained using the formula [4} 
These Q factors, as a function of OD are 
shown in figures [4] and [5| 

The above data was all fit to a to a 
stretched exponential function 

Q = A ex p((x — XqY + c, 


. Here, A is the amplitude, x is the Optical 
Density, x 0 is the shift in scaling and /3 is 
the stretch coefficient. The dependency is 
normal exponential for smoke, whereas it is 
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Figure 4: Variation of Q factors with respect to OD 
for smoke, vertical/horizontal and 45/135 degrees 



Figure 5: Variation of Q factors with respect to OD 
for fog, vertical/horizontal and 45/135 degrees 


a stretched exponential with stretch coeffi¬ 
cient of f3 ~ 0.4 for vertical and horizon¬ 
tal polarizations, that too in case of smoke. 
The variation is not surprising since smoke 
and fog particles have different distributions 
in terms of size and scattering cross sec¬ 
tions. Since the OD is a quantitative mea¬ 
sure of number of scattering events per unit 
length, a higher OD corresponds to more 
scattering and therefore a stronger depolar¬ 
ization. As seen from the graphs, Q fac¬ 


tor decreases as OD increases and asymp¬ 
totically reaches zero. However, The exact 
physics behind this behaviour is needed to 
be explored further. Part of the answer may 
lie in the calculations of Moeyaert et. al. 
ra, who have shown that the BER in case 
of an OOK has a stretched exponential de¬ 
pendency. 

In case of 45°/ 135°, the variation of Q 
factor with respect to OD is slower com¬ 
pared to the ‘H/V’ basis, with both po¬ 
larization modes showing a stretched expo¬ 
nential dependency with Q = Aexp[—(x — 
xo) 0 ' 5 , where x is the optical density with 
a shift of x 0 . This is very critical in 
Quaternary encoding, since the 45°and 135 
° polarized are depolarized more than the 
vertical/horizontal polarizations. However, 
it is evident that the differential method of 
computing SOP and then mapping to infor¬ 
mation bits has a higher tolerance for depo¬ 
larization. 

Bit Error Rate(BER) 

Theoretical estimate of BER, based on Q 
factor is given by equation [4| which is more 
valid for a OOK scheme However in 

this case we can define two BER values, one 
for each polarization - 

The estimated BER values are very close 
to zero for OD almost upto -20 dB, after 
which it raises sharply. The actual val¬ 
ues of BER, which are obtained by compar¬ 
ing transmitted and received bits are shown 
in [Tj These indicate a very similar trend, 
though the error rates for smoke are much 
lower than that for fog. This is due to the 
fact that fog causes a significantly higher 
depolarization than smoke. 

These results can be used to estimate a 
threshold for noise tolerance in using real 
communication situation. Instead of mea¬ 
suring actual bit error rates in real time, 






























(a) 



(b) 



Figure 6: Estimated Bit Error rates in presence of 
smoke. Solid line indicates an exponential fit with 
stretch factor. Values are reasonably low for almost 
upto OD=-20 dB. 



OD OD 



(a) 




OD OD 

(b) 

Figure 8: BER for 45° (left) and 135° (right), as a 
function of signal attenuation, in presence of (a) fog 
and (b) Smoke. Solid line is fit for exponential func¬ 
tion. As in case of vertical/horizontal case, BER is 
close to zero till a threshold value of OD, and in¬ 
creases sharply beyond that. But this increase is 
much sharper and higher for smoke increases much 
more than in case of fog. 


4. Stokes Vector analysis 


Figure 7: Experimentally obtained BER from com¬ 
paring transmitted and received data in presence of 
fog (left) and smoke (right). Solid lines are fit for 
exponential curve. 

one can rely on the signal attenuation and 
propose that if the signal strength is less 
than -15 dB, the communication channel is 
likely to be noisy. An attenuation of more 
than 20 dB will definitely be unreliable. 

It was shown in reference m that the 
value of SOP is determined only by snake 
and ballistic photons which retain their po¬ 
larization and the depolarized part of the 
light will cancel out and does not affect the 
SOP. We show below the same results using 
a full Stokes vector analysis. 


Stokes vector picture gives a complete 
analysis for this situation since the depolar¬ 
ization due to atmosphere is properly repre¬ 
sented by the corresponding Muller matrix. 
Stokes vector then gives a correct represen¬ 
tation of a partially polarized light. 

We obtain the Muller matrices for all the 
optical components used, as M\ for the first 
PBS, M 2 for the depolarizing field in free 
space and the M 3 for the second PBS. Since 
PBS has two outputs, the corresponding 
Muller matrices are represented by Mj and 
M/, with the primed components standing 
for horizontal and unprimed components in¬ 
dicating outputs at vertical port of PBS. 
They are given as 


6 









































































( 6 ) 
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0 

2 
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0 

0 
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and 


M[ = M 3 



( 1 

-1 

0 

°\ 

1 

-1 

1 

0 

0 

2 

0 

0 

0 

0 


^0 

0 

0 

0/ 


where, a is the depolarization factor. For 
an ideal case of no depolarization, a = 1 
and a = 0 for a complete depolarization. 
In this situation, the outputs of polarizing 
beam splitter (M 3 and M 3 ) are equivalent to 
taking projections of the input Stokes vector 
onto two different states, each of which gives 
a probability of 0.5. 

Net measurement from detectors Di and 
D 2 would then be 


D 1 = (M 3 .M 2 .M{)L 2 
D 2 = (M! i .M 2 .M[)L 2 


Depending upon whether LI or L2 is 
pulsed. Photodiodes measure total inten¬ 
sity of light, which corresponds to So of the 
Stokes Vector. It can be shown from from 
equations [5] and [ 6 j that 
for LI case 




1 + Cl\ 
1 T a 
0 


D 2 = 


\ 0 / 

and for L2 case 


1 - a \ 

— 1 — CL 

0 

0 / 


(7) 


D\ = 


1 — a \ 
— 1 — a 
0 

0 ) 


Do = 


1 d* a \ 
— 1 — a 
0 

0 ) 


( 8 ) 


Using this in[TJ for only the So component 
D 1 = (M 3 .M 2 .M 1 )L 1 + (M 3 .M 2 .M[)L 2 of the Stokes vector ’ we § et for L1 case as 


D 2 = (M' z .M 2 .Mi)Li + (M^.M 2 .M[)L 2 

Li = {1,1,0,0} and L 2 = {1,-1,0,0} 
are the Stokes vectors for light from lasers 
LI and L2. During the communication ei¬ 
ther LI is fired or L2 is fired. The two 
lasers are never simultaneously on. Hence 
the stokes vector of the light at the detec¬ 
tors reduces to 


D\ — (iff 3 . M 2 .M\) Li 

D 2 = (M^.M 2 .Mi)Li (5) 

or 


{Di}i — {D 2 } 1 _l + a — l + a_ 

{Di}i + {D 2 }i 1 + a + 1 — a 

and the LI case as 

{D 2 } 1 — {Di}i _ 1 — a — 1 — a_ 

{D 2 }i + {Di}i 1 — cl -j-1 + a 

( 10 ) 

It is evident from above that the value of 
SOP is either —a or + 0 . With no depolar¬ 
ization this will be ± 1 , and reduces to lower 
values when a depolarization occurs. If the 
light is completely depolarized, then a = 0 
and SOP can not be determined. 


7 




















Similarly, for the case of 45°and 135 “case, 
we can add the Muller matrix for half wave 
plate as well 


1 0 0 0\ 

0 0 10 

0-100 

0 0 0 1 / 

The Stokes vector on the detector would 
therefore be 



D\ — (^M^.M 2 .M 4 

D 2 = (M^.M 2 .M 4 )L 1 (11) 

D\ = i y M^.M 2 .M 4 )L 2j 
and 

D 2 = -\-(^M^.M 2 .M 4 )L 2 (12) 


respectively. The Corresponding Stokes 
vectors, now in 45°and 135°polarized light 
would be Li = {1,0,1,0} and the L 2 = 
{1,0, —1, 0}. It can be shown that the final 
Stokes vector for light falling on the detec- 


tors would be 




/l + a\ 


/ 1 — a \ 


1 + a 

0 

d 2 = - 
2 

— 1 + a 

0 


1 0 ) 


1 0 / 


(13) 


This gives an SOP as 
{A}i — {D 2 }i _ 1 + a — 1 + a 


{fdi} 1 + {-A}! 1 + a + 1 — a 

Similarly when laser L 2 is on, we get 
{A}i — {-A}i —1 + a + 1 + a 


= a (14) 


{Ah + (Ah 


— 1 + a — 1 — a 


= —a 


( 15 ) 


As in case of vertical polarized fight, the 
SOP values are either —a or +a, which can 
be mapped to digits zero and one. This 
mapping is reliable as long as a is not equal 
to zero. 

5. Conclusion 

We have proposed a Quaternary encoding 
using four states of polarization, for use of 
free space Optical communication. To un¬ 
derstand the nature of depolarization due to 
atmospheric effects such as fog and smoke, 
we created these effects within laboratory 
situation and studied it. It shows that be¬ 
haviour of depolarization is different for fog 
and smoke, due to the difference in their 
constituent nature, but also for vertical and 
45°basis. This can cause serious errors in a 
communication protocol. 

However, the Quality factors and Bit er¬ 
ror rates obtained from a statistical distri¬ 
bution of the states of polarization, show 
a reasonably high Quality factor and a low 
BER, even for a significantly high scattering 
densities. In addition, a differential mea¬ 
surement method to map the final bit to the 
polarization allows even lower noise rates. 
The results shown will be very important in 
a practical implementation of a free space 
Quaternary encoding communication setup, 
which will be our plan for further studies. 
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